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1 Introduction 

1> . The parameter \Vub\ is one of the smallest and least known elements of the Cabibbo- 



^ ! Kobayashi-Maskawa (CKM) quark- mixing matrix. A precise determination of \Vub\ 

O I would significantly improve the constraints on the unitarity triangle and provide a 

^ ■ stringent test of the Standard Model mechanism for CP violation. With the CKM 

■ angle 03, \Vub\ can constrain the unitarity triangle from tree level processes alone. 

^ . Experimental studies of charmless semileptonic B decays can be broadly catego- 

^ ! rized into inclusive and exclusive measurements depending on how the final states 

I I are treated. The inclusive method measures the decay rate r(i? — > Xuii'), where Xu 

is known as the hadronic system that does not contain charm-quark. On the other 
1^ • hand, the exclusive method measures the decay rates for exclusive final states such 

a.s B ^ Tfiu and piu. Two methods give not only different efficiencies and signal-to- 
^ ! background ratios, but also different theoretical calculations to be used in order to 

^ I extract \Vub\- Using both approaches and comparing the results will help us verify 

the robustness of the theoretical error estimation, which dominates the current un- 
certainty in the determination of \Vub\- Progress in last few years will be summarized 
in this presentation. 



2 Inclusive determination of \Vub\ 

The inclusive semileptonic decay rates in the quark level depend only on CKM matrix 
element and the quark mass, as shown in the following equation 

T{b^uiV) = -^^\V^b\'ml (1) 
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where is the Fermi constant and rrib is the 6-quark mass. The hadronic level is 
easy to calculate with the framework of operator product expansion (OPE) [1] and 
the \Vub\ can be parametrized as 

|K.|= 0.00424| ^(^^^^"^'^^^-^^^P^ |^ X (1.0±0.012qcd±0.022hqe) (2) 

where the first error comes from the uncertainty in the calculation of quantum chro- 
modynamics (QCD), perturbative and non-perturbative quantities and the 2nd from 
the uncertainty in the heavy quark expansion (HQE), sensitive to nib. This formula- 
tion |2j has been updated in ICHEP06. The main problem in the inclusive method 
is the background from the 6 to c decays because the rate is approximately 50 times 
larger, namely 

T{b-^ciu) ^50T{b-^uiiy). (3) 

One may attempt to remove this large background by applying kinematic selection 
criteria but once the signal to background ratio becomes controllable, the OPE is 
known to fail in such limited phase space. In the theory side, in order to overcome 
this problem, various different techniques have been developed. For example, a non- 
perturbative shape function [3] (SF) is developed to extrapolate to the full phase 
space. The shape function is the lightcone momentum distribution function of the 
6-quark inside the meson. The detailed shape is not known theoretically from the first 
principle and even phenominologically, the low-tail part is least known. The shape 
function is needed to be determined from experimental data. The other approach is 
called dressed gluon exponentiation (DGE) [1] according to the factorization proper- 
ties of the fully differential width in inclusive decays Sudakov logarithms exponentiate 
in moment space. The third approach [5] measures the ratio of |Kife|/|VJs| with the 
photon energy spectrum in 6 57 decay mode. This technique is interesting as the 
only residual shape function dependence remain at the end. 

• Lepton Endpoint Analysis 

The theoretical calculations allow for the extraction of the observed partial B X^ii^ 
decay rate above a certain lepton momentum to the total inclusive B Xud-v decay 
rate using the measured shape function parameters and a subsequent translation of 
the total decay rate to \Vuh\- The experimental study was pioneered by the CLEO 
collaboration ^ and the recent work was carried out by the BaBar collaboration 
Fig.[T](a), (b), and (c) show the electron energy spectra for various cases. Open circles 
in Fig. [T] (a) show the on-resonance data and closed circles with a curve in Fig. [T] (a) 
show the off-resonance data where non-i?i? background is included. The triangles in 
Fig. [1] (b) show the data after the non BB background subtraction and the histogram 
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in Fig. [T] (b) show the simulated BB background. Closed squares and the histogram 
in Fig. [1] (c) show the data and the simulated signal after all background subtraction, 
respectively. It is obvious that the subtraction of backgrounds is extremely crucial in 
this kind of analysis. The shaded region in Fig. [1] (c) is used for the final extraction 
of the \Vub\- In the mean time, the partial branching fraction can be obtained from 
the background subtracted data and the summary from three experiment is listed in 
Table [TJ Note that the Belle's result has the lowest cut on Ei. 




Electron Momentum (GeV/c) 



Figure 1: The distribution of electron momentum is shown |7j. Open circles in 
(a) show the on-resonance data and closed circles with a curve in (a) show the off- 
resonance data where non-BB background is included. The triangles in (b) show 
the data after the non BB background subtraction and the histogram (b) show the 
simulated BB background. Closed squares and the histogram in (c) show the data 
and the simulated signal after all background subtraction, respectively. 



• Measurements of mx, P+, and 

In this analysis, the measurements are made with a sample of events where the 
hadronic decay mode of the tagging side B meson, Btag, is fully reconstructed, while 
the semileptonic decay of the signal side B meson, Bgig, is identified by the presence 
of a high momentum electron or muon. B denotes both charged and neutral B 
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Table 1: Summary of the partial branching fraction oi B ^ Xuijy given the lepton 
energy cut. 



data 


Ai3(10~^) 


El (GeV) 


CLEO 9/fb 


2.30 ± 0.15 ± 0.35 


2.1 


Belle 27/fb 


8.47 ± 0.37 ± 1.53 


1.9 


BaBar 80/fb 


5.72 ± 0.41 ± 0.65 


2.0 



mesons. This method allows the construction of the invariant masses of the hadronic 
(Mx) and leptonic system in the semileptonic decay, and the variable P+ = 

Ex — \px\ where Ex is the energy and the magnitude of the three-momentum 
of the hadronic system. These inclusive kinematic variables can be used to separate 
the B — ^ X^iu decays from the much more abundant B — >■ Xcii' decays. Three 
competing kinematic regions were proposed by theoretical studies [31 [8], based on the 
three kinematic variables, and are directly compared by this analysis. 
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Figure 2: The P+ distribution (symbols with error bars) after subtracting B X^iiy 
background, with fitted B X^fz/ contribution (histogram). 

The value of \ Vub\ is extracted using recent theoretical calculation [3l[8] that include 
all the current known contributions. Results from the Belle experiment using a data 
set of 253/fb [9J is the first one to use P+. Figure [2] show the distribution of the 
variable Pj^. Signal region is defined by Px < 0.66 GeV and shows a clear indication 
of enhancement of the signal in Fig. [2l BaBar did also similar analysis but only with 
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Table 2: Summary of the partial branching fraction oi B ^ Xuijy given the lepton 
energy cut. 



data 


Phase Space 


AB (10-4) 




Mx < 1.7 


12.4 ± 1.1 ± 1.0 


Belle 253/fb 


Mx < 1.7, g2 > 8 


8.4 ± 0.8 ± 1.0 




P+ < 0.66 


11.0 ± 1.0 ± 1.6 


BaBar 211/fb 


Mx < 1.7, g2 > 8 


8.7 ± 0.9 ± 0.9 
(preliminary) 



Mx and [TO]. The measured partial branching fractions from both experiments 
for different phase space values are summarized in Table [2l Note that the errors are 
larger than those appeared in the endpoint analyses results. 

• Extraction of \Vub\ 

In order to extract the \Vub\ from the partial branching fraction measurements 
described so far, one has to know the distribution of the shape function. The photon 
energy spectrum in i? — > Xg'j provides access to such distribution function of the 
b quark inside the B meson [11]. The knowledge of this shape function is a crucial 
input to the extraction of \Vub\ from inclusive semileptonic B X„£z/ measurements. 
Both Belle and Babar fit the spectrum to theoretical predictions in order to extract 
the \Vub\. 

Using the heavy quark parameter values m6(SF) = 4.60 ± 0.04 GeV and /i^(SF) 
= 0.20 ± 0.04 GeV^, the value of \Vub\ from various experimental results is extracted 
within "BLNP" framework [3] . This is done by the heavy flavour averaging group [12] 
(HFAG) and the result is 

IK^IbLNP = (4-49 ± 0.19exp ± 0.27theory) >< 10^'- (4) 
Note that the total error in percentage is 

^iKfelBLNP = ±7.3% (5) 

and the result indicates that the precision of the measurement is still limited by 
theoretical uncertainty but the value of the error is not far from the experimental 
error. The individual result from experiments is shown in Fig. [3l Note that the 
best measurement is from the endpoint analysis at present. Contributions to the 
experimental and theoretical errors can be found in Table [31 The contribution from 
sub-leading SF, for example, is known to be hard to reduce from the present value. 
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CLEO (endpoint) 
4.09 ± 0.48 ± 0.36 

BELLE (endpoint) 
4.82 ± 0.45 ± 0.30 

BABAR (endpoint) 
4.39 ±0.25 ±0.39 
BABAR (E„ q-) 
4.57 ±0.31 ±0.41 
BELLE 

4.06 ± 0.27 ± 0.24 
BELLE sim. ann. (mj,, q^) 
4.37 ± 0.46 ± 0.29 
BABAR (mx, q') 
4.75 ±0.35 ±0.32 
Average +/- exp +/- (mb, theory) 
4.49 ±0.19 ±0.27 

X"/dof=6.1/6(CL = 40.7%) " 

OPE-HQET-SCET (BLNP) 

Phys.Rev.D72;073006,2005 

in,j input from b-4 c 1 V and b— > s y moments 



HFAG 



Figure 3: Summary of the measurements of \Vub\ within BLNP framework from 
various experiments are shown. The averaged value is iKiblgLNp = (4.49± 0.19exp± 

0-27theory) xlO"'- 



One can also extract \Vub\ using the DGE framework [4J. This is also done by 
HFAG and the result is 

\Vub\BGE = (4-46 ± 0.20exp ± 0.20theory) >< 10"'- (6) 

with the parameter mf,(MS) = 4.20 ± 0.04 GeV. Note that the result is in good 
agreement with the result with BLNP method, and is remarkable as two theoretical 
methods use rather different approach in their calculations. The contribution of the 
error is also listed in Table HI The sharing of the error is also similar to the case of 
the BLNP framework. One thing to note is that the weak annihilation is not taken 
into account in this case. 

On the other hand, BaBar explored alternative methods in obtaining \Vub\- Lei- 
bovich. Low, and Rothstein (LLR) [5] have presented a prescription to extract \Vub\ 
with reduced model dependence from either the lepton energy or the hadronic mass 
mx The calculations of LLR are accurate up to corrections of order and 
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Table 3: Summary of the contributions to the experimental and theoretical errors of 
\Vui,\ in the BLNP framework. 



Source 


contribution (%) 


statistical 


2.2 


Expt. systematic 


2.8 


b — > civ model 


1.9 


b uiu model 


1.6 


HQ parameters 


4.2 


sub-leading SF 


3.8 


Weak Annihilation 


1.9 



Table 4: Summary of the contributions to the experimental and theoretical errors of 
|T4b| in the DGE framework. 



SoTirce contribution (%) 



statistical 


1.8 


Expt. systematic 


2.5 


b — > cliy model 


2.3 


b — > uiiy model 


2.3 


rub (R_CUT) 


1.2 


as (R_CUT) 


1.0 


Total semiloptonic width 


3.0 


DGE theory 


2.9 
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[KniB / {Cf^b)Y 1 where is the experimental maximum hadronic mass up to which 
the B Xud-v decay rate is determined and A ~ ^QCD' This method combines 
the hadronic mass spectrum, integrated below with the high-energy end of the 
measured differential B Xg'j photon energy spectrum via the calculations of LLR. 
The measured | Vub I as a function of ( is shown in Fig. |H The small arrow indicates 
the value of ( that is used for the cut, ( = 0.67 GeV/c^. At this point, the measured 
value becomes 

IK^I = (4.43 ± 0.38stat ± 0.25^j,t ± 0.29^1,) x IQ-^ (7) 




1 2 3 

^ [GeV/c^] 

Figure 4: \Vub\ as a function of ( with the LLR method (left) and for the determi- 
nation with the full rate measurement (right). The error bars indicate the statistical 
uncertainty. They are correlated between the points and get larger for larger ( due 
to larger background from B Xclv. The total shaded are illustrates the theoret- 
ical uncertainty; the inner light shaded area indicates the perturbative share of the 
uncertainty. The arrow indicates C, = 1.67 GeV/c^. 

Another approach that BaBar chooses to reduce the model dependence is to mea- 
sure the B —>■ Xuii' rate over the entire Mx spectrum. Since no extrapolation is 
necessary to obtain the full rate, systematic uncertainties from and Fermi mo- 
tion are much reduced. Perturbative corrections are known to order a^. The rate 
of -B ^ Xuii' is extracted from the hadronic mass spectrum up to C = 2.5 GeV/c^ 
which corresponds to about 96 % of the simulated hadronic mass spectrum, and find 
\Vub\ = (3.84 ± 0.70stat ± 0.30syst ± O.lO^j^) x 10"^ using the average B lifetime 
of Tb = (1.604 ± 0.012) ps. The current uncertainties on the B Xg'-f photon en- 
ergy spectrum limit the sensitivity with which the behavior at high ( can be probed. 
These two new results are consistent with previous measurements but have substan- 
tially smaller uncertainties from rrib and the modeling for Fermi motion of the b quark 
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inside the B meson. Both techniques are based on theoretical calculations that are 
distinct from other calculations normally employed to extract \Vub\ and, thus, provide 
a complementary determination of \Vub\- 



3 Exclusive determination of \Vub\ 

In this section, we discuss the extraction of \Vub\ from exclusive decays such as B 
{tt, p,u!)ii'. For B^ —> 7i~£~^i' decays, the differential decay rate becomes 



247r3 



where f+iq"^) is a form factor and is the squared invariant mass of the i'^u system. 
Only shape of f+{q^) can be measured experimentally. Its normalization is provided 
by theoretical calculations which currently suffer from relatively large uncertainties 
and, often, do not agree with each other. As a result, the normalization of the 
f+{q^) form factor is the largest source of uncertainly in the extraction of \Vub\ from 



the B^ 



TT 



V branching fraction. Values of f+iq"^) for B^ 



TT 



v decays are 



provided by unquenched [HI [15] and quenched [T6j lattice QCD calculations, presently 
reliable only at high q^ (> 16 GeV^/c"^), and by Light Cone Sum Rules calculations [H] 
(LCSR) based on approximations only valid at low q^ (< 16 GeV^/c^), as well as by 
a quark model [18]. The QCD theoretical predictions are at present more precise for 
B^ -K^l^v decays than for other exclusive B X^lv decays. Experimental data 
can be used to discriminate between the various calculations by measuring the f+{q^) 
shape precisely, thereby leading to a smaller theoretical uncertainty on \Vub\- 




Figure 5: Fit to dV /dq^ on exclusive B ^ n £+z/ decays. 
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For the recoiling the other side of the B meson, we may require nothing (untagged) 
or require indication of semileptonic decay (D^*Hu tag), and require hadronic decay 
(full reconstruction tag). We review recent progress on experimental studies with 
different tagging methods listed above. 

• Untagged B niu 

CLEO pioneered the measurement of Vub with exclusive decays |l19J . They perform 
a simultaneous maximum likelihood fit in AE and Mmiu to seven sub-modes: tt^, tc^, 
p^, uj/r] pi^n^Ti^ , and f] — > 77. In the fit they used isospin symmetry to constrain 
the semileptonic widths r^^^^i) ^ 2T^^ (vr^) and r^^^^i) ^ 2T^^ (p°) ~ 2T^^ 
(u), where the final approximate equality is inspired by constituent quark symmetry. 
Signals for vr and p are extracted separately in three bins. Given form factors from 
theory, they extracted \Vub\ from a fit to dT/dq'^, and it is shown in Fig. [51 From 
their fit, we see that ISGW2 [18] is least favored. Combining B nil' and B piu 
results, CLEO found that 



iKfel = (3.17 ± 0.17(stat)l:|]:}^(syst)lgj^(theo) ± 0.03pp) x 10^1 (9) 

BaBar also did similar but improved analysis [20]. More accurate value of 
was obtained in the so-called F-average frame where the pseudo-particle Y has a 
four- momentum defined by Py = {Pn + Pe)- The angle 6yt between the directions 
of the p*^ and py momenta in the T{4S) rest frame can be determined assuming 
energy-momentum conservation in a semileptonic B Yu decay. The use of the 
F-averaged frame yields a q^ resolution that is approximately 20 % bettern than 
what is obatined in the usual T(4S') frame where the B meson is assumed to be at 
rest. They also fit the /S.B/B spectrum using a probability density function based on 
the f+{q^-, o) parametrization of Becirevic-Kaidalov [2T] (BK). The normalized /S.B/B 
distribution is shown in Fig. [6l together with the result of a f+i^f) shape fit using 
the BK parametrization and theoretical prediction. The obtain a value of a = 0.53 
± 0.05 ± 0.04. BaBar data are clearly incompatible with the ISGW2 quark model, 
which is in agreement with what CLEO data indicated before. The extraction of \Vub\ 
is carried out from the partial branching fractions using \ Vub\ = A.B/{r^AQ, where 

= (1.536 ± 0.014) ps is the B^ lifetime and A( is the normalized partial decay rate 
predicted by various form factor calculations. For the LCSR calculations with < 
16 GeVVc^ \Vub\ = (3.6 ± 0.1 ± 0.1 1^) x 10"=^ is obtained. For the HPQCD and 
and FNAL lattice calculations with q^ > 16 GeVVc^ |Kfe| = (4.1 ± 0.2 ± 0.2 to.l) 
X 10~3 and |Kfe| = (3.6 ± 0.2 ± 0.2 l[Jj) x 10"=^ are obtained, respectively. Note 
that all are in good agreement within given uncertainty and this gives us confidence 
in exclusive measurement of \Vub\- 
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Figure 6: Differential decay rate formula fitted to the normalized partial ISBjB 
spectrum in 12 bins of g^. The smaller error bars are statistical only while the larger 
error bars include statistical and systematic uncertainties. The BK parametrization 
(solid black curve) reproduces the data quite well (x^=8.8 for 11 degrees of freedom) 
with the parameter a = 0.53 ± 0.05 ± 0.04. The data are also compared to LCSR 
calculations (dotted line), unquenched LQCD calculations (long dashed line and short 
dashed line) and the ISGW2 quark model (dash-dot line). 



• L)(*)£z/ tag 

Belle presented measurements of 5° ti~ /p~l^v and n^/p^i'^v decays 

using B — ^ D^*^^^ tagging p2|. They reconstruct the entire decay chain from the 
T(4S') -^sig-^tag' -^sig ~^ 7r/p£z/ and E^^^g ~^ D'^*Hi' tag with several D^*) sub- 
modes. The back-to-back correlation of the two B mesons in the T(4S') rest frame 
allows us to constrain the kinematics of the double semileptonic decay. The signal 
is reconstructed in four modes, B^ — > -n^ jp^V^v and B^ -n^jp^t^v. Yields and 
branching fractions are extracted from a simultaneous fit of the B^ and B^ samples 
in three intervals of g^, accounting for cross- feed between modes as well as other back- 
grounds. Belle applied this methods to 5 — >• n j piv decays for the first time, and have 
succeeded in reconstructing these decays with significantly improved signal-to-noise 
ratios compared to the //-reconstruction method. With the data of 253 fb~^. Belle 
extracted branching fractions and |\4fe|. Figure [7] shows the distribution for the 
decay B^ — > it' /p~t^v. The error bars are too large to reject any of form factor 
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Table 5: Summary of \Vub\ measurements by the Belle collaboration with different 
lattice calculations. 





g2 GeVVc^ 






FNAL 


> 16 


3.60 ± 0.41 ± 0.20 ± 


0.62 
-0.41 


HPQCD 


> 16 


4.03 ± 0.46 ± 0.22 ± 


0.59 
-0.41 



models for the moment. Table [5] summarized the results with two different lattice 
calculations. This gives about 13 % experimental uncertainty on \ Vub\, currently dom- 
inated by the statistical error of 11 %. By accumulating more integrated luminosity, 
a measurement with errors below 10 % is feasible. With improvements to unquenched 
LQCD calculations, the present method may provide a precise determination of \Vub\- 
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Figure 7: Extracted distribution for — TT^i'^u decays. Data points are shown 
for different form factor models used to estimate the detection efficiency. Lines are 
for the best fit of the form factor shapes to the obtained distribution. 



The BaBar collaboration also did similar analysis |23] based on 211 fb~^ of data. 
They also included the hadronic decay of B mesons for tagging. The detailed tech- 
niques are similar to what to be explained later. The summary of the measurement of 
\Vub\ can be found in Table El All the results are in good agreement with the results 
from the Belle experiment. 

• Full reconstruction tag 

In this method [25], Belle fully reconstructs one of the two B mesons from T(4S') 
decay (-Stag) following hadronic decay modes, B^ D'^*^^Tr^, B^ 
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Table 6: Summary of \Vub\ measurements by the BaBar collaboration with different 
theoretical calculations. 





g2 GeVVc^ 


\Vub\x 10-3 




Ball-Zwicky 


< 16 


3.2 ± 0.2 ± 0.1 


+U.5 
-0.4 


HPQCD 


> 16 


4.5 ± 0.5 ± 0.3 


+0.7 
-0.5 


FNAL 


> 16 


4.0 ± 0.5 ± 0.3 


+0.7 
-0.5 


APE [23] 


> 16 


4.1 ± 0.5 ± 0.3 


+1.6 
-0.7 



D^*^^ai, or B^ Ag -^^as done in other B meson analyses, decays are 

identified on the basis of the proximity of the beam-energy constrained mass M-^^^ 
and AE to their nominal values of the B meson rest mass and zero, respectively. If 
multiple tag candidates are found, the one with values of M-^^ and AE closest to 
nominal is chosen. Events with a -Stag satisfying the selections M-^^ > 5.27 GeV/c^ 
and —0.08 < AE < 0.06 GeV are retained. The charge of the -B^ag candidate is 
necessarily restricted to Qtag = or Q^ag 1 by demanding that it is consistent with 
one of the above decay modes. Reconstructed charged tracks and electromagnetic 
clusters which are not associated with the -Stag candidate are used to search for 
the signal B meson decays of interest recoiling against the -Btag- After all cuts to 
enhance the signal are applied, remaining data are projected as a function of missing 
mass squared (^niiss)' shown in Fig. [HI As is demonstrated in the figure, a high 
purity signal can be obtained with this method. A preliminary branching fractions 
are obtained as 

B{B = 1.49 ± 0.26(stat) ± 0.06(syst) x 10-^ (10) 

B{B ir^iu) = 0.86 ± 0.17(stat) ± 0.06(syst) x 10"^ (11) 

Whilst the statistical precision of these measurements is limited at present, the po- 
tential power of the full reconstruction tagging method, when it can be used with 
larger accumulated 5-factory data samples in the future, can clearly be seen. 

4 Conclusions 

We have discussed recent progress in measurements of \Vub\- At present, the total 
error from the inclusive measurement is approximately 7 %. The theoretical uncer- 
tainty in it is still larger than the experimental uncertainty, but not by a lot anymore. 
The shape function, describing the Fermi motion of the 6-quark inside the meson still 
remains a big issue in extracting \Vub\, but different approaches produce consistent 
numerical values indicating the problem is well understood. Also, BaBar's new ap- 
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Figure 8: Missing mass squared (^^jgg) distributions after all cuts, for B — > n^Eu 
on the left and B n^ii' on the right. Data is indicated by the points with error 
bars. The blue histogram (lightest shade in greyscale) show the fitted prediction 
based on the LCSR model. The green histogram (middle shade in greyscale) shows 
the fitted b uiu background contribution. The crimson histogram (darkest shade 
in greyscale) shows the fitted background contribution from other sources. 



proach based on Leibovich, Low and Rothstein [5] may look promising in future as it 
has residual SF dependence only. 

The exclusive measurements of \Vub\ gives the total error to be greater than 10 % 
at present and the form factor is the main issue in this field. The untagged analyses 
are the most precise ones at present. The theoretical uncertainty is still dominant in 
exclusive measurements. 

I am grateful to Masahiro Morii from the BaBar experiment who kindly provided 
complete information on their experimental results. 
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